Adiabatic Cooper-Pair Pumping in a Linear Array of Cooper Pair Boxes 
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We present a study of adiabatic Cooper pair pumping in one dimensional array of Cooper pair 
boxes. We do a detailed theoretical analysis of an experimentally realizable stabilized charge pump- 
ing scheme in a linear array of Cooper pair boxes. Our system is subjected to synchronized flux and 
voltage fields and travel along a loop which encloses their critical ground state of the system in the 
flux-voltage plane. The locking potential in the sine-Gordon model slides and changes its minimum 
which yields the Cooper pair pumping. Our analytical methods are the Berry phase analysis and 
Abelian bosonization studies. 
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1. Introduction: The adiabatic pumping physics 
gets more attraction after the pioneering work of Thou- 
less [J, Q ■ Quantum adiabatic pumping physics is related 
with many systems like open quantum dots [H, 0, [Bj] , su- 
perconducting quantum wires [fl LZlK J osephson junctions 
[il>0] > the Luttinger quantum wire |lfj |. interacting quan- 
tum wire 



Uj and also to the quantum spin pump 



quar 
0. 



An adiabatic parametric quantum pumping is a device 
that generates a dc current by a cyclic variation of sys- 
tem parameters, the variation being slow enough that 
the system remains close to the ground state throughout 
the pumping cycle. It is well known to us that when a 
quantum mechanical system evolves then it acquires a 
time dependent dynamical phase and time independent 
geometrical phase [l3| . Geometrical phase dependent on 
the geometry of the path turned in the parameter space. 
For the closed path encircling the critical ground state 
differ by the phase factor and the cyclic parametric vari- 
ation introduce non-vanishing transport in the system. 
Quantum pumps are to transform an ac signal of fre- 
quency / into dc current given by the relation / = fQ, 
where Q — ne for electron pumping and Q = 2ne for 
Cooper pair pumping. Quantum pumping with perfect 
accuracy could be utilized to establish a standard of cur- 
rent. Here we have considered the perfect pumping con- 
dition. The error in the pumping procedure arises due 
to the current reversal and the spontaneous charge exci- 
tations. In Ref. ([3] ), the authors have discussed the 
sources of errors (the non-adiabatic correction leave the 
system in an unknown superposition of the charge state, 
instead of definite charge state, after the full cycles) and 
their minimization. Here we would like to explain the 
theoretical detailed of that experimental proposals stabi- 
lized charge pumping 14[ in a linear array of Cooper pair 
boxes. Before we start our full swing quantum held theo- 
retical calculations. We would like to present a derivation 
to illustrate the Cooper pair pumping in one dimensional 
array of Cooper pair boxes. 



2. MODEL HAMILTONIANS AND 
CONTINUUM FIELD THEORETICAL STUDY: 

Our Hamiltonian (H = H c + Hj) consists of two 
parts one is the Josephson coupling and the other 
is the Coulomb charging energy. Where Hj = 



EfcLi E Jik cos(f) k and H c = J2k=i let' ® k is the cnar S e 
occurs at the kth junction, this charge is measured w.r.t 
the gate charge. So the effective charge of the k'th junc- 
tion can be varied by varying the gate voltage of the 
system connected to the junctions. The charging Hamil- 
tonian is diagonal in the basis formed by the charge eigen 

state |n >= \n\,n2, n>N-i >> n % is the number of 

Cooper pair in each island. We would like to explain the 
transport of Cooper pair of this one dimensional tunnel 
junction through the analysis of Berry phase [13j. We will 
see that the Cooper pair current in the system consists 
of two parts one is the conventional supercurrent and the 
other is adiabatic Cooper pair pumping. Here we follow 
the seminal paper of Berry [f3| for our theoretical analy- 
sis. In our system, we are varying the Josephson junction 
couplings through the applied flux and the charge on the 
dot by applying the gate voltage. So the adiabatic vary- 
ing parameters are, R = (Ej, Qk). The state \ip(t) >, of 
the systems evolves according to Schrodinger's equation 

H(R(t))\ip{t) >=ih\ip(t) > . 

At any instant, the natural basis consists of the eigen- 

states \n(R) >= \ni R ,n 2 R, n N -x R >. \tp(t ) >= 

\n(t ) >. 
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m) > = < 

< l(t )\V R n(t ) > .S t R\l(t ) > 
where F(E u E n ) = ( e - iE ^ - e~ iE ^) 



(1) 



\iP(t) >= \n(t) > + \Sn{t) > 



(2) 
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FIG. 1: The stabilized Cooper pair pumping system. Identical 
pulse sequence is applied to every alternate superconducting 
island. SET is the single electron transistors, they are biased 
with short voltage pulses every time the pump should be in 
a definite charge state. Leftmost and rightmost island are 
coupled to the reservoirs, m is the number of Cooper pair at 
the ith site. V g and C g are respectively the gate voltage and 
capacitance of each Cooper pair box. Circle with two crosses 
represents the Josephson junction. 



coupling and the electrostatic potential of superconduct- 
ing island in a Cooper pair The adiabatic 
quantum Cooper pair pumping procedure, is nothing but 
the transport of Cooper pair from one end of the system 
to the other end. The basic of Cooper pair tunneling in 
an array of Cooper pair boxes can be understood from 
the analysis of Cooper pair transport in a single Cooper 
pair box with two terminal Squid. One of the terminal 
(say left) to transport a Cooper pair into the box and 
another to transport the same pair to the reservoir. This 
process generate the current. The electrostatic energy of 
the system can be expressed as 



i 



lg(i+l), 



(8) 



V_r is the derivative w.r.t externally varied parameter. 
The amount of charge that passes through the tunnel 
junction k during St interval of time is 



SQ k 



< tl>(t)\Is,k\il>(t) > dt 



(3) 



SQ k = 5t < I S , K >„(t ) +2Re[ / < n(t)\I s , k \Sn{t) >]. 

Jto 

(4) 

Where 



Here E c u\ and E m ^ ii+ i^ are respectively the charging 
energy and the electrostatic couplings between two is- 
lands. We are interested in the charge degeneracy point, 
i.e., when the gate charge is close to 1/2, the lowest en- 
ergy states are characterized by either zero or one Cooper 
pair on each island. With this assumption they have 
reduced the Hilbert space and map the system to a fi- 
nite anisotropic Heiscnberg spin-1/2 chain in an external 
magnetic field In the spin language Cooper pah- 

pumping is nothing but the transport of spin (Jordan- 
Wigner fermions) from one end of chain to the other end. 
They have defined the Hamiltonian [3] 



N 



Is,k = -£'52Ej ik (-i\n R + 6 k ><n R \e i *+h.di) (5) H = -±B x l o x x -\b x * 1 a x N ' - £ 



Sk measure the change of n due to the tunneling of one 
Cooper pair through the kth junction. The first term 
gives the direct Cooper pair transform, i.e., the super- 
current . The second term corresponds to the adiabatic 
Cooper pair transform. The charge transfer contribution 
for the second term is 



SQi = 



-2hTlJ <nlIs ^ l> < llSn> ] (6) 



l=£n 



El — E n 



The total charge transport due to Cooper pair transport 
over a pumping period r is given 



Q 



-2e 
~7T 



d 4> E n (t)dt+2hY,Im[— n ^ S ' k 



>< l\Sn > 



Ei — E n 



(7) 

We are interested to build up the theory of experimen- 
tal proposal of Ref. [14J. There system consists of N 
Cooper pair boxes. Nearest neighbor are coupled with 
dc- SQUIDS. Leftmost and rightmost island are coupled 
to the reservoirs via other Squids as shown in Fig. 1. This 
figure is nothing but the collection of tunable Josephson 



N-l 

+- [A. M+ ia z V/ +1 - J M+1 (a+V_ J+1 + a+ i+1 aj§] 
i=i 

Ai^+i is the constant electrostatic coupling amplitude. 
The tunable parameters of the system: B X 1,N the Joseph- 
son coupling of the leftmost and rightmost Squids. B z l is 
the electrostatic potential of the island and Ji^+i is the 
Josephson coupling between the neighboring island. We 
are interested in the charge degeneracy point at this point 
the most favorable state of the system is the antiferro- 

magnetic configuration ( 1 010101 > and 1 101010 > 

). They start with one of the antiferromagnetic states 
and transfer the charge of every island to the right by 
two sites to achieve pumping. They have implemented 
it by applying identical pulse sequence to every second 
island [14]. So we can write B x x = Ja^ = = J e , Q and 
B X N — J\fl = J3 7 4 = J 0je . Je j0 and J 0ie are respectively 
the even-odd and odd-even Josephson couplings. They 
have also considered a difference between the charging 
energies between the odd and even sites. In our theo- 
retical analysis, we consider the Josephson couplings for 
even and odd sites as respectively Ej 1 = Ej (l — Si(t)) 



3 



and Ej 1 — Ej (1 + 8\ it) ) . The charging energies of even 
and odd sites are respectively B z \ — Sq(1 — <Sa(i)) and 
B z2 = B (l + (t) ) ■ One can write the Hamiltonian in 
terms of spin operators: 



H = -E^o(i-(-ir<M*))(s+"s- 



n+l 



S+ n+1 S- 



n n 

One can express spin chain systems to a spinless fermion 
systems through the application of Jordan- Wigner trans- 
formation. In Jordan- Wigner transformation the rela- 
tion between the spin and the electron creation and 
annihilation operators are S% = ipfyn ~~ 1/2 , S~ — 
ip n exp[iirY^Z- 00 n j] > S i = e M~i^T,"=-oo n ]} > 
[l8| . where rij = ipltpj is the fermion number at site j. 



n 

+A^(^,i - l/2)(^ +1 V'n+i - 1/2) , 

-f 1 - {-iTHtmUn 1/2). (ii) 

n 

Our approach is completely analytical. Before we pro- 
ceed further for continuum field theoretical study of these 
model Hamiltonians, we would like to explain the basic 
aspects of quantum Cooper pair pumping in terms of 
spin pumping physics of our model Hamiltonians: An 
adiabatic sliding motion of one dimensional potential, 
in gapped Fermi surface (insulating state), pumps an 
integer numbers of particle per cycle. In our case the 
transport of Jordan- Wigner fermions (spinless fermions) 
is nothing but the transport of spin from one end of 
the chain to the other end because the number opera- 
tor of spinless fermions is related with the z-component 
of spin density We shall see that non-zero S(t) 

introduces the gap at around the Fermi point and the 
system is in the insulating state (Peierls insulator). In 
this phase spinless fermions form the bonding orbital 
between the neighboring sites, which yields a valance 
band in the momentum space. It is well known that the 
physical behavior of the system is identical at these two 
Fermi points. One can analyses this double degeneracy 
point, from the seminal paper of Berry [l3j |. It appears 
as source and sink vector fields defined in the generalized 
crystal momentum space. B n (K) = V k x A n (K), and 
A n (K) = ± < n(K)\V K \n(K) >, where K = (M(t)). 
Here B n and A n are the fictitious magnetic field (flux) 
and vector potential of the nth Bloch band respectively. 
The degenerate points behave as a magnetic monopole 
in the generalized momentum space fl3| , whose magnetic 
unit can be shown to be 1, analytically 13, IH 



dS ■ B± = ±1 



(12) 



positive and negative signs of the above equations are re- 
spectively for the conduction and valance band meet at 
the degeneracy points. Si represent an arbitrary closed 
surface which enclose the degeneracy point. In the adia- 
batic process the parameter S(t) is changed along a loop 
(r) enclosing the origin (minima of the system). We de- 
fine the expression for spin current (I) from the analysis 
of Berry phase. Then according to the ori gina l idea of 
quantum adiabatic particle transport 0,11 IH, E3, the 



total number of spinless fermions (J) which are trans- 
ported from one side of this system to the other is equal 
to the total flux of the valance band, which penetrates 
the 2D closed sphere (S2) spanned by the T and Brillioun 
zone [12j. 



1=1 dS ■ B+i = 1 
's 2 



(13) 



We have already understood that quantized spinless 
fermion transport is equivalent to the spin transport . 
We will interpret this equation more physically after Eq. 
18 . This quantization is topologically protected against 
the other perturbation as long as the gap along the loop 
Studies of spin pumping explain 
pumping against 



una 



remains finite 

the stabilization of quantized spin 
z-component of exchange interactions. 



We recast the spinless fermions operators in terms of 
field operators by this relation 



i>{x) 



_ I e ik F x 



xl> R (x) + ^ lkFX ^l(x)} (14) 



where ipn( x ) an d describe the second-quantized 

fields of right- and left-moving fermions respectively. We 
want to express the fermionic fields in terms of bosonic 
field by this relation 



-i (rcf>(x) - 9(x)) 



V2t 



(15) 



r is denoting the chirality of the fermionic fields, right (1) 
or left movers (-1). The operators U r are operators that 
commute with the bosonic field. U r of different species 
commute and U r of the same species anticommute. <f> 
field corresponds to the quantum fluctuations (bosonic) 
of spin and 9 is the dual field of (p. They are related by 
this relation 4>r = — <f> and 4>l = + 4>- 

U sing the standard machinery of continuum field the- 
ory [18( , we finally get the bosonized Hamiltonians as 



Ho = v 



[ L ^{^ : n 2 
Jo 2^ 



[dx^x)] 2 :} 



2A 



4A 



dx : {d x Mx)f : + : [d x <t> R {x)? ■ 



si 



+ —j- / dx {d x (f) L {x)){d x (j) R (x)) 



(16) 
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Hq is the gapless Tomonoga-Luttinger liquid part of the 
Hamiltonian with vo = sinkp- 

After continuum field-theory the Hamiltonian become 



P s z increases by 1 per cycle. We define it analytically as 



H 



H, 



EjJijt) 
2n 2 a 2 
B 6 2 {t) 
2ixa 



dx : cos[2VK(/>(x)) 



dx : cos[2v K<fi(x)] 



- 9 9 / dx : cos[aVK(/)(x)} : - — 
lir z a z J I 



dxdr. 



The second term of the Hamiltonian for NN exchange 
interaction has originated from the XY interaction. This 
dimerization is the spontaneous dimerization, i.e., the in- 
finitesimal variation of Josephson coupling in lattice sites, 
is sufficient to produce a gap around the Fermi points. 
The third term of the Hamiltonian arises due to the site 
dependent on-site charging energies modulated by the 
gate voltage. It yields the staggered phase of the sys- 
tem. The effect of applied gate voltage on the Cooper 
pair box appears as an effective magnetic field and also 
as a staggered magnetic in the spin representation of the 
model Hamiltonian. The system is in the mixed phase 
when both interactions (2nd and 3rd terms of the Hamil- 
tonian) are in equal magnitude otherwise the system is 
in any one of the states of the mixed phase depending 
on the strength of the couplings. The last term can be 
absorbed in the Hamiltonian through the proper shifting 
of the wave function. So when 1/2 < K < 1 only these 
(2nd and 3rd terms of the Hamiltonian) time dependent 
dimerizing field (third term of Eq. 17) is relevant and 
lock the phase operator at <fi = + . Now the locking 
potential slides adiabatically (here the cyclic magnetic 
flux and gate voltage fields that produces the dimeriza- 
tion). Speed of the sliding potential is low enough such 
that system stays in the same valley, i.e., there is no scope 
to jump onto the other valley. The system will acquire 
2-7T phase during one complete cycle of dimerizing field. 
This is the basic mechanism of Cooper pair (spin pump- 
ing) of our system. The quantized Cooper pair transport 
of this scenario can be generalized up to the value of A for 
which K is greater than 1/2 . In this limit, z-component 
of exchange interaction has no effect on the Cooper pair 
pumping physics of Hamiltonian. This expection is easily 
verified when we notice the physical meaning of the phase 
operator (<jf> (x)). Since the spatial derivative of the phase 
operator corresponds to the z-componcnt of spin density 
(Cooper pair density) , this phase operator is nothing but 
the minus of the spatial polarization of the z-component 
of spin, i.e., P s z = — jj J2jLi J^j Z - Shindou has shown 
explicitly the equivalence between these two considera- 
tion [l2j]. During the adiabatic process < (fit > changes 
monotonically and acquires - 2tt phase. In this process 



5P S 



dP s Z = / d.ri), 



x) >= 1 



(18) 



This physics always hold as far as the system is locked by 
the sliding potential and A < 1 [l2| ■ The above equation 
for Cooper pair transport is physically consistent with 
the Eq. 13 (based on Berry phase analysis) of Cooper 
pair current. The quantized Cooper pair transport of 
this scenario can be generalized up to the value of A for 
which K is greater than 1/2 . In this limit, z-component 
of exchange interaction has no effect on the Cooper pair 
pumping physics of our system. 

In summary, we have presented the theoretical explana- 
tion of adiabatic Cooper pair pumping of experimentally 
reliazable stabilized charge pumping scheme of an array 
of Cooper pair boxes. The charge state of our system is 
definite. 

The author (SS) would like to acknowledge The Center 
for Condensed Matter Theory of IISc and also Mr. Va- 
sudeva for reading the manuscript very critically. 
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